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Abstract. We study the dispersion relations and spectra of in- 
variant Schrodinger operators on a graphyne structure (lithographite) . 
In particular, description of different parts of the spectrum, band- 
gap structure, and Dirac points are provided. 



Introduction 

Graphene, a monolayer of graphite, is famous for its unusual elec- 



tric and mechanical properties (e.g., 10,19]). Recently, researchers 
suggested other 2D carbon allotropes which were given the common 
name "graphynes." It has been suggested (see, e.g. (2|[T8]) that some 
graphynes, which have not been synthesized yet, might be even more 
interesting than the graphene. 

Various standard and less standard approaches have been used to 
model the spectral structure of graphene and graphynes (one of the 
most popular was to use a version of density functional technique [Tj). 
One of the ways similar 2D structures have been modeled previously, 
was using the techniques of quantum networks, also known as quan- 
tum graphs (see, e.g., [l||2^). In particular, several studies of spectra 
of Schrodinger operators on graphene and carbon nanotube structures 
(e.g., 11,12,17]) have been conducted, which have proven to be much 



simpler to study and preserving all essential ingredients of the disper- 
sion relation. One should also be aware of a recent study of disper- 
sion relations of 2D Schrodinger operators with honeycomb symmetry 
in [8|[9], where in particular the mandatory presence of Dirac cones is 
established. 

In this paper, we take the quantum graph approach similar to 17 to 
study spectra of Schrodinger operators on the simplest graphyne among 
14 various structures suggested in ff] (it represents the 2D projection 
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of the SO called litliographite [5]). From now on we reserve the word 
"graphyne" for this particular structure. We derive the dispersion re- 
lations for these operators on graphyne. From here, we extract various 
information about the spectral structure of the operators. Unlike sim- 
ilar periodic operators in M", the quantum graph operators can (and 
often do) have point spectrum (i.e., bound states). We find this part of 
the spectrum and provide an explicit description of the corresponding 
eigenspace^ The presence of spectral gaps and conical "Dirac" points 
is also studied. The formulations of the results involve the discriminant 
of the Hill operator with the potential obtained by periodic extension 
of the ID potential on a single edge. 

In Section [T] we introduce the geometry of the structure and the 
operators of interest. In Section |2] we derive the dispersion relation 
and the band-gap structure for graphyne with the main results stated 
in Theorems [7] and [9] The proof of an auxiliary Proposition 10 is given 
in Section m 



1. Geometry of graphyne structure and related 
schrodinger operators 

The graphyne structure that we study is represented by the graph 
G shown in Fig. [T] and has square symmetry, unlike the graphene's 
honeycomb one. At each vertex there is a carbon atom that is bonded 
to other atoms. Chemical bonds between atoms are represented by the 
edges connecting the corresponding vertices. 



All the edges of G are assumed to have length 1. We denote by E{G) 
and V{G) correspondingly the set of all edges and all vertices of G. 
There is a free action of the group 1? of integer vectors in M? on G 
by the shifts by vectors piCi + ^262, where (^1,^2) ^ and ei = 
(\/3, 0),e2 = (0,2). We choose the domain W shown in Fig. [l]as the 
fundamental domain of this action. It contains three vertices a, 6, c 
and pieces of five edges /, h, k, I as shown in the figure. We choose 
the directions of these edges as shown. Notice that the boundary of 
the chosen fundamental domain does not contain any vertices (which is 
always possible to achieve). The reason for using such domain is that 
the presence of a vertex on its boundary would unnecessarily complicate 
considerations. The entire structure G can be obtained from W by Z^- 
shifts, which also define directions on all edges of the graph G. 



The presence of bound states is an artefact of the quasi-lD model. It, however, 
often indicates possible presence of very flat bands in the "grown up" system. 
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Figure 1 . The structure G and a fundamental domain 
W with vertices a, 6, c and parts of edges /, g, h, k, I (and 
some of their lattice shifts) 

We define on each directed edge e the arc length coordinate Xe that 
identifies it with the directed segment [0,1]. When it does not lead 
to ambiguity, we will use x instead of Xe to denote the coordinate on 
the edge e. One can introduce now in a natural way the Hilbert space 
L2(e) as the space of all square integrable functions on the edge e and 
H2{e) as the Sobolev space on e that consists of functions with two 
distributional derivatives in L2(e). 

We also define 

L^iG) = L,{e) 

eeE(G) 

as the space of all square integrable functions on G. 

We will use the notations Ue for the restriction of a function u on G 
to an edge e. We also use u'^ for the derivative of Ue in the direction of 
the edge e. 

Let qo{x) be an even and real L2-function on [0,1], i.e. qo{x) = 
go(l — x) for a.e. x G [0, 1]. 

Using the described before identification of the directed edges with 
the segment [0,1], we can transfer the potential qo{x) to each edge, 
thus defining a potential q{x) on the whole G. 

It is not hard to show that the evenness assumption on go implies 
the following property: 

Lemma 1. The potential q defined as above is invariant with respect 
to the symmetry group of the graph G. 

We are now ready to construct the Schrodinger operator H in L2{G), 
whose spectral properties will be studied in this paper. The operator 
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Hu{x) = 3—^ h q{x)u{x). 



H acts on each edge as follows: 

d'^u{x) 
dx"^ 

Its domain D{H) consists of all functions u{x) on G such that: 
(1) Me G H2{e), for all e € E(G'), 
(2) 

Yl \MH2ie) < (2) 

eeE(G) 

(3) at each vertex these functions satisfy Neumann vertex condition, 
i.e. Mei(f) = Me2('^^) ^OT any edges 61,62 containing the vertex v 
and 

= 0, for any vertex v in V{G). 

Thus defined operator is well-known to be unbounded and self-adjoint 
(e.g. [3| Theorem 1.4.19]). It is also invariant with respect to all sym- 
metries of the graph G. 

2. Graphyne spectrum 

In this section we study the spectrum of the operator H. Let us 
describe first the main steps of our approach. The technique of Floquet- 
Bloch theory [3||6|[T5|[2l] allows us to reduce the consideration to a 
family of spectral problems on the fundamental domain W . Then one 



can switch to a discrete problem (e.g., 20 , [3| S ection 3.6]). This step 
uses the standard Hill's operator theory [6||21|. Finally, the discrete 
problem can be analyzed rather explicitly. 

Let us get to some detail now. For each Q = {Q\, 62) in the Brillouin 
zone B = [— 7r,7r]^, let be the Block Hamiltonian that acts as ([l]) 
on the domain that consists of functions u{x) that belong to Hf^J^G) 
and satisfy Neumann vertex condition along with the following cyclic 
(or Floquet) condition: 

u{x + P161 + ^262) = u{x)e'^'^ = u{x)e'^P'^'+P^'^''\ (3) 

for all {pi,P2) G and all x G G. 

Due to this condition, such functions u are uniquely determined by 
their restrictions to the fundamental domain W. 

We have the direct integral expansion § Section 4.3] 

H= H^de. 

Jb 
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Therefore jsj Theorem 4.3.1], 

a{H)= U a{H'). 

6»e[-7r,7r]2 

It is well-known (e.g., |[3|,fT3|) that operator has purely discrete 
spectrum (t{H^) = {Aj(6')} with lim Xj{6) = oo. The multiple- valued 

j-s>oo 

function 6 i — > {Xj{(^)} is called dispersion relation and its graph - 
dispersion surface or Block variety of the operator H. Spectrum of H 
is the range of the dispersion relation for 6 changing in the Brillouin 
zone. Thus, we now concentrate on studying the spectrum of H^, for 
6 E B, i.e on solving the eigenvalue problem: 

H% = Am, A G M, (4) 

for u G H2iW) satisfying the cyclic condition ^ at the boundary and 
Neumann vertex condition inside. 

Combining vertex and cyclic conditions we have: 

' Uf{0) =Ug{l) =Uh{0) =Uk{l) =: A 

n'/0)-n;,(l) + <(0)-<(l) = 

Ufil)=uiil) = uh{l)e''^ =:B 
' w'/l)+izKl)+<(l)e^'^=0 

%(0) = UkiO)e''' = uiiO)e-''^ =: C 
^ u'g{0) + <(0)e*^i + u'i{0)e-'^^ = 0. 

We will need another auxiliary operator: 

Definition 2. We denote by the Dirichlet Hamiltonian on [0, 1] 

that acts as ^ with Dirichlet boundary conditions u(0) = u{l) = 0. 
We also denote by the (discrete) spectrum of . 

For each A ^ S^, there exist two linearly independent solutions 
(po,(pi such that ipo{0) = ipi{l) = l,ipo{l) = <^i(0) = 0. Sometimes we 
address (po, ipi as ipo^x and ipi^\ to emphasize their dependence on A. 
We use the same notation ipQ, ipi for analogous functions on each edge 
of W under fixed identification of these edges with the segment [0, 1], 
which should not lead to a confusion. Then for A ^ solution of (H) 
can be represented as follows: 

Uf = AifQ + Bifi 
Ug = CipQ + Aipi 
< Uk = Ce~''^^(fo + Aifi 
Uh = Aipo + 5e-^^Vi 
ui = Ce^^Vo + Bifi. 
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Continuity and eigenvalue equation on each edge are already satisfied. 
What is left to be checked is the zero flux condition at each of the three 
vertices in W: 







M(o)(i 



35^;(l) + Ce^^Vo(l) =0 
5e-*^V;(0) + 3C/o(0) = 0. 



(5) 



Notice that ip[{l) = -</?o(0) and ip[{0) 
of function go- Thus ([s]) becomes 



-(/3q(1) due to the evenness 



+ B^[{0){1 + e-''^) + C(^;(0)(1 + e- 
Aip[{0){l + e*^i) - 3Bip[{l) + Ce*^2<^;(0) = 
A(/?;(0)(1 + e*^i) + 5e-^^V'i(0) - 3^/^(1) = 0. 



(6) 



Since <f'i{0) 7^ 0, we can define 



<a(0)' 



(7) 



then dGI) is reduced to 



-4r/(A)A + (1 + e-'^')B + (1 + e-'^')C = 
(1 + e'^')A - 3r7(A)S + e^'^^C = 
(1 + e'^')A + e~'^^B - 3r/(A)C = 0. 

Determinant of this system is 

-4[97]^i\) - r]{\) - (cos^i + l)(3r7(A) + cos^s)]- 
These calculations prove the following: 

Lemma 3. A point A ^ is in the spectrum of the Schrddinger 
operator H if and only if there exists 9 = {61,62) G B such that x = 
ri{X) is a root of the equation 

- X - (cos ^1 + 1) (3x + cos ^2) = 0. (8) 

Let us now extend go periodically from [0, 1] to the whole real axis 
M and consider the Hill operator H'^'^'" on M as below: 

i7P-«(x) = -^^ + go(x)n(a;). 

Here we use the same notation go(a;) for the periodic extension. The 
monodromy matrix M(A) of H^^^ is defined by the formula 



M(A) 



>(0)" 
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where satisfies the differential equation 

- + qo{x)^{x) = X^{x) on M. (9) 

Discriminant (or Lyapunov function) trM[X) of the Hill operator H^^^ 
is denoted by D{\). Next proposition ( [T7| Proposition 3.4]) collects 
some well-known results about the spectra of Hill operators \6\: 

Lemma 4. 

(1) The spectrum alH^'^^) of H^'^^ is purely absolutely continuous. 

(2) aiRP"') = {A e M||D(A)| < 2}. 

(3) a{HP^^) consists of the union of closed non- overlapping (al- 
though, possibly touching) and non-zero length finite intervals 
(bands) ■= [^2^, &2fc], ^2^+1 := [&2fc+i, ^2^+1] such that 

ao < bo < hi < ai < a2 < b2 < ■ . ■ 

and lim = 00. 

k—^oc 

The (possibly empty) segments (62^, ^2fc+i) (ind {a2k,a2k+i) O'l"^ 
called the spectral gaps. 

Here {0^} and {bk\ are the spectra of the operators with periodic 
and anti-periodic conditions on [0, 1] correspondingly. 

(4) Let G be the /c*^ Dirichlet eigenvalue labeled in increasing 
order. Then belongs to (the closure of) the k*^ gap. When 
Qo is even, \^ coincides with an edge of the k-th gap. 

(5) If X is in the interior of the k*^ band B^, then -D'(A) 7^ 0, 
and -D(A) is a homeomorphism of the band onto [—2,2]. 
Moreover, -D(A) is decreasing on (— oo,&o) CL'^d {a2k,b2k) and is 
increasing on (&2fc+i; '22A:+i)- It has a simple extremum in each 
spectral gap [0^,0^+1] and [bk,bk+i]- 

(6) The dispersion relation for H^'^^ is given by 

D{X) = 2cos^, 

where 9 is the one- dimensional quasimomentum. 

Claim (4) of the lemma about the even potential case can be ex- 
plained as follows: let u{x) be the eigenfunction of the Hill operator 
ffper corresponding to the fc*'* Dirichlet eigenvalue Af", i.e. 

HP^'u{x) = X^u{x),u{0) = u{l) = 0. 

Then u{l — x) is also an eigenfunction corresponding to that eigenvalue. 
Either u{x)-\-u{l — x) or u{x)—u{l — x) is nonzero and therefore will be 
an eigenfunction corresponding to A^. Since u{x) — x) is periodic 
and u{x) — u{l — x) is anti-periodic, A^ must coincide with an edge of 
the k^'^ gap. 
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The following relation between function ri{X) (see ([7])) and the dis- 
criminant -D(A) of i/P^'" is well-known (6|[T7] and easy to establish: 

V{X) = l-D{X). (10) 



Therefore, according to the statement (2) of Lemma |4] and equality 
(10), one needs to analyze the roots of the cubic equation (jsj). 

So far we have just dealt with A not in the Dirichlet spectrum of 
. Now let us consider the exceptional values A G S^. 

We introduce the following notion first: 

Definition 5. An eigenfunction is said to be a simple loop state if it 
is supported on a single hexagon or rhombus of the structure G and 
vanishes at all vertices (see Fig. 

We can now describe what happens for A G S^. 

Lemma 6. Each A G T,^ is an eigenvalue of infinite multiplicity of the 
operator H . The corresponding eigenspace is generated the simple 
loop states. 

Proof. For each A G S"^, let ipx be the corresponding eigenfunction 
of operator . Since go is even, one can assume function ipx to be 
either even or odd. For an odd eigenfunction ipx., we repeat it on each 
edge of hexagon/rhombus; for even eigenfunction ipx - repeat around 
hexagon/rhombus with an alternating sign. In both cases we get an 
eigenfunction of H that lives only on one particular loop. Thus A G 
app{H). Let us prove infinite multiplicity of the eigenvalues A G S^, 



Figure 2. Simple loop states constructed from even 
function on [0, 1] for hexagon and odd function on [0, 1] 
for rhombus 

which is a well-known feature of periodic problems. Let Mx C L2{G) 
be the corresponding eigenspace and 7 be some period vector of G. The 



^I.e., is the closed linear hall of 
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shift operator by 7 acts in Mx as an unitary operator. Suppose that 
Mx is finite dimensional, then has an eigenfunction / G Mx C L2{G) 
corresponding to an eigenvalue fi such that = 1. On the other hand, 
/ is multiplied by /i when shifted by the vector 7. Since = 1, / 
clearly cannot belong to L2{G), which leads to contradiction. Thus 
A G are eigenvalues of infinite multiplicity. 

In order to prove that the eigenspace is generated by simple loop 
states with hexagonal and rhomboidal supports, it is enough to prove 
that these simple loop states generate all compactly supported eigen- 
functions in the eigenspace Mx- Indeed, as it is shown in [m] (see 
also [sj Theorem 4.5.2]), linear combinations of compactly supported 
eigenfunctions are dense in the space Mx- 

First we notice that each compactly supported eigenfunction of 
H vanishes at all vertices. Indeed, due to connectedness of G, there 
must be a "boundary" vertex v of the support that is connected by an 
edge with a vertex w outside the support. We claim that {p{v) = 0. 
Otherwise, we have an edge such that the function vanishes at one end, 
w (corresponding to x = 0) and does not vanish at the other end. We 
introduce a basis of solutions of (|4]), functions cx and sx, such that 
ca(0) = 1,sa(0) = sa(1) = (a non-trivial function sx exists, since 
A G T,^). The eigenfunction ip can be represented as ip{x) = Acx{x) + 
Bsx{x). In particular, = ip{0) = A, and so 7^ (/?(!) = Bsx{l) = 0, 
which leads to contradiction. Repeating this argument, we conclude 
that the eigenfunction ip vanishes at all vertices. Besides, the support 
of if cannot have a vertex of degree 1. (Otherwise, due to Neumann 
boundary condition, both function and its derivative will vanish at that 
vertex, which makes function to be equal to zero.) 

Now one needs to prove that (p can be represented as a combination of 
simple loop states. Consider the external boundary of the support of ip, 
which is a closed circuit C of edges, containing the whole support inside. 
The interior of this curve is a union of N elementary hexagons and/or 
rhombuses of the graph G. We begin with a boundary edge Cq G G. 
One of the N internal hexagonal or rhomboidal loops must contain cq. 
Let ipo be the simple loop state that coincides with (p on the edge Cq and 
is extended to that loop as described before. Function ip — ipo will be 
the new eigenfunction with a smaller support (number of loops — 1). 
Continuing this process (see Figjs]), we will eventually represent the 
eigenfunction (p as a combination of simple loop states. □ 

In the next theorem, which is our main result of this section, we 
describe the dispersion relation and the structure of the spectrum of 
operator H. 
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Figure 3. An example of deleting simple loop states 
(the dark ones) from the support of an eigenfunction 

Let F{6) be the triple-valued function providing for each 6 the three 
roots of the equation ([s]). By Proposition [Io| which we will formulate 
and prove later, function F is real- valued in the Brillouin zone. Assume 
F{9) = {Fi{9),F2{e),F3{e)), where Fi{9) < FiiO) < F^iO) for all 9 E 
B. Then we have 

Theorem 7. 

(1) The singular continuous spectrum asc{H) is empty. 

(2) The dispersion relation of operator H consists of the following 
two parts: 

i) pairs {9,\) such that 0.5D{\) e F{d) (or, A G D-\2F{9))), 

where 9 changing in the Brillouin zone; 

and 

a) the collection of fiat (i.e., 9 -independent) branches A G T,^ . 

(3) The absolutely continuous spectrum aac{H) has band-gap struc- 
ture and is (as the set) the same as the spectrum a{HP'^^) of the 
Hill operator H^'^^ with potential obtained by extending periodi- 
cally go from [0, 1]. In particular, 

aUH) = {XeR\\D{X)\<2}, 

where D{X) is the discriminant of H^'^^ . 

(4) The bands of cx^H) do not overlap (but can touch). Each band 
of a{H^'^^) consists of three touching bands ofa{H). 

(5) The pure point spectrum (Tpp{H) coincides with T,^ and belongs 
to the union of the edges of spectral gaps of a{H^^'^) = aac{H). 
Eigenvalues A G S-^ of the pure point spectrum are of infinite 
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multiplicity and the corresponding eigenspaces are generated by 
simple loop (hexagon or rhombus) states. 
(6) Spectrum (j{H) has gaps if and only if a{H^'^'^) has gaps. 

The statements of the theorem are illustrated in Fig. |4j 



bands gaps 




Figure 4. The bold segments are the bands of a{HP^'^). 
Each of them is split into three touching bands of cr(if). 
One eigenvalue at the end of a band is also shown. 



Proof. The first claim about emptiness of the singular continuous spec- 
trum of Schrodinger operator H is rather well-knowirj 

Let first A ^ S^. Then, according to Lemma pT {9, A) is in the 
dispersion surface of H iff 77(A) is a root of ([s]) for this 9. In other 
words, due to (10), 0.5-D(A) must be one of the three values of F{9). 
Hence, this is equivalent to D{X) e 2F{9), or A G D-\2F{9)). 

If A G S^, then according to Lemma [6| (^, A) is in the dispersion 
surface for any 9 from the Brillouin zone. 

This proves the second statement of the theorem. 

According to the Lemmas |4] and |6] we have 

C o-(/7),S^ C (T(ifP"") (11) 



and 

For A ^ 
for some 9. 
of equation 

"'it goes back to the famous L. Thomas' absolute continuity theorem l24| 
e.g., Isl Theorem 4.4.1], or 21 and references therein. 



^^Jjper^ = {A G M| \D{X)\ < 2}. 

A is in the spectrum of H iff //(A) is a root of equation ^ 
Proposition 10 below shows, in particular, that all roots 
^ belong to [—1, 1] and cover this interval. Thus, A is in 



See, 
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the spectrum of H iff \rj{\)\ < 1. Since -D(A) = 2?7(A), this means that 
a{H) \ = a{HP^'') \ and, by closure, a{H) = a^H^^''). 



The same Proposition 10 shows that the graph of the triple-valued 
function F{9) does not have any flat branches outside the set of values 
S^. Thus, the spectrum of is absolutely continuous outside S^. This 
argument, together with Lemma [6] finishes the proof of the statements 
(3) and (5) of the theorem. 

From statement (2), the dispersion relation of H consists of the va- 
riety A = D~^{2Fj{6)),j G 1,3, and collection of fiat branches A G 
located at some edges of spectral bands. 

According to Lemma |4| function -D(A) is a monotonic homeomor- 
phism from each spectral band of the Hill operator onto [—2,2]. Be- 
sides, the ranges of functions 2Fi, 2F2 and 2F3 are all in [—2, 2]. Thus, 
part of the spectrum a{HP'^^) that corresponds to each band of Hill 
operator ffP^'" coincides as a set with the part of the absolutely contin- 
uous spectrum aac{H) that consists of three bands. In another words, 
operator H has "three times more" non-fiat bands than Hill operator 
HP"'' does. 

The function D^^ is multiple- valued, thus producing infinitely many 
bands for any 2Fj,j = 1,3. Two such spectral bands (pre-images 
of the same 2Fj{B)) clearly cannot overlap, because function -D(A) is 
monotonic on each band. Two spectral bands which are pre-images of 
2Fj{B) and 2Fi{B) for different i,j G 1,3 also cannot overlap because 
the ranges of functions Fj,j = 1,3 belong to [—1, 1] and do not overlap 



by Proposition [10] below. One should notice that although spectral 
bands do not overlap, they still can touch each other, and as we will 
see below this is indeed the case. This can happen at points {9, A) 
such that -D(A) = ±2 or ±2/3. This proves the statement (4) of the 
theorem. 

Since the spectra of H and H^^^ coincide as sets, we get the last 
statement of the theorem. □ 

Corollary 8. 

(1) Unless the potential go is constant, the spectrum cx^H) has at 
least one gap. 

(2) For a generic smooth potential qo, all possible gaps in cr{H) are 
open. 

Proof. The first claim of the Corollary follows from the last statement of 
Theorem [7] and Borg's theorem 13. Similarly, the second claim follows 
from Simon's genericity result l22] instead of Borg's theorem. □ 
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Graphene has captured physicists' interest because of its unusual 
electronic properties. These properties are caused by occurrence of so- 
called conical singularities or Dirac points. Roughly speaking, Dirac 
points are points where two spectral bands touch and locally form a 
cone (also known as Dirac cone). One is interested in conical singu- 
larities that are stable under small perturbation of the potential not 
breaking the symmetry. 

In the next theorem, we will take a closer look at the spectral bands 
of the operator H. Moreover, we will specify all the conical singularities 
in the Brillouin zone B and describe how the spectral bands behave 
near these points. 

In what follows, we will use the notation 

^0 arccos(— 1/3). 

Theorem 9. 

(1) In the free case, i.e., when the potential qq is equal to zero, 
the Bloch variety of H has conical singularities at the following 
points: 

i) {e, A) = (0, 0, {2{k + l)7r)2), at which D{X) = 2, 

ii) [e, A) = (0, ±7r, ((2A; + l)7rf), at which D{X) = -2, 
Hi) {e, A) = (±^0, 0, (^0 + 2A;7r)2), at which D{X) = -2/3, 

iv) {e, A) = (±^0, ±7r, (27r - + '^knY), at which D{X) = 2/3, 
fork = 0,1,2,... 

(2) When we turn on a small potential qo ^ 0, the conical singular- 
ities corresponding to |-D(A)| = 2 will generically split into two 
smooth branches, opening a gap. The conical singularities that 
occur when |-D(A)| = 2/3 are stable under small perturbation by 
a potential q of the type considered. 

Before proving this theorem, we need to explore some properties of 
function F{9). 

Proposition 10. 

(1) Function F(9) is real-valued for 9 in the Brillouin zone B. 

(2) Let F{e) = {F,{9),F2{9),F3{9)), where F^{9) < F2{9) < F,{9) 
for all 9 E B. Then the ranges of functions Fi,F2,Fz are 
[—1, —1/3], [—1/3, 1/3] and [1/3,1] correspondingly. 

(3) Function Fi attains its maximal value at {9i, 0) for9i G [— tt, — ^o]U 
[^o-tt] or (±7r, ^2) for ^2 G [— tt, tt] and minimal value at 9 = 
(0,±7r). 

Function F2 attains its maximal value at (^i,±7r) for 9i e 
[—^0) ^0] aiT'd minimal value at 9 = {9i, 0) for 9i e [—9o, 9o\. 
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bands gaps 




Dirac points 



Figure 5. The Dirac points are located at the points 
where three bands of (t{H) touch to form a band of 

Function F3 attains its maximal value at (0, 0) and minimal 
value at (^i,±7r) for 61 G [— vr, — ^0] U [^O;^] or (±71,6*2) for 
62 e [-7r,7r]. 

(4) The linear level sets of function F{6) inside B are 
{(±7r,^2),^2e [-7r,7r]}, 
{(^i,±7r),^iG[-7r,7r]}, 
{(ei,±7r/2),^^i G [-7r,7r]}, 
{(0i,O),^iG[-7r,7r]}. 

Function F{6) does not have any flat branches. 

Proof of this proposition will be given in Section [3] 
Proof of Theorem 

1. In the free case, D{X) = 2cosV^ (see, e.g., IgI), so we have 



cos^/X = Fj{e),j = 1,2,3. (12) 

As it was proven in Theorem [7| spectral bands do not overlap. It is 
still possible that these bands touch each other at their edges. We will 
study now whether this indeed happens and prove that at these points 
the spectral bands have conical form. 



According to (12), all non-flat spectral bands are 

Agfc+j = (arccos(Fj(6')) + 2A;7r)^, Agfc+s+j = (27r — arccos(Fj(0)) + 2A;7r)^, 

for all j = l^,k = 0,1,2,. . . 
Thus for = 0, 1, 2, . . ., we have 

i) Bands Xek+4, Aefc+r touch each other at (0, 0, (2(fc+l)7r)^), for which 
D{X) = 2. 

ii) Bands X^k+s, ^ek+e touch each other at (0, ±7r, ((2A; + l)vr)^), for 
which D{X) = -2. 
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iii) Bands \ek+2, Aefc+a touch each other at (±6*0, 0, {9o + 2A;7r)^) while 
bands X^k+n, Aefc+e touch each other at (±6*0, 0, (— ^^o + {2k + 2)Tr)^). 
At these points D{X) = -2/3. 

iv) Bands X^k+i, XQk+2 touch each other at (±6*0, ivr, (vr — 6*0 + 2A;7r)^) 
while A6A:+4, Aefc+s touch each other at (±^o, ±7r, {71 + 60 + ^kn)^). 
At these points ^(A) = 2/3. 

Let us now look at the structure near the touching points. One needs 
to deal with each case above separately. Since the argument we use for 
i) and iii) are similar to those needed for ii) and iv), for the sake of 
brevity, we will consider only cases i) and iii). 

i) Let {9, A) = (0, 0, Aq) where Aq = {2{k + l)7r)2, keN. Then 



D{X) 



cos 



2 

for A — 7- Aq where a2 = — I/SAq < 



VX = 1 + a2(A - Ao)2 + o((A - Aq)'), (13) 



and 



q2 

cos^i = 1 - ^ + o{el), for 01^0 (14) 



02 

cos^a = 1 - Y + o{dl), for 02 0. (15) 



Since 0.5D{X) is a root of the equation ([8|, we have 



D{X)Y D{X) . ^ ,.(D{X) 



9^-^ ^ = (cosei + l) 3^-^ + cos^2 • (16) 



Plugging expressions ( |13[ ), ( |14[ ) and (15) into (16) and simplify the 
expression, one obtains the following formula: 

A{X - Xof + 261 + dl = o{dl) + o{dl) + o((A - Ao)'), (17) 



for (^^, A) — (0,0, Ao) where A = 20a2 < 0. Equation (17) shows that 
the spectral bands touching at the point {9, A) = (0, 0, Aq) have the 
conical form. In other words, the Bloch variety of operator H has 
conical singularity at (6*, A) = (0, 0, Aq). 

iii) Let now {6, A) = (^o, 0, Aq) where Aq = {eo + 2k7c f, A; = 0, 1, 2, . . .. 
We have 

^ = + ai(A - Ao) + a2(A - Ao)' + o((A - Ao)'), 
for A — i- Ao, where ai = 0.5D'{Xo), 

cos^i = -l + biiOi - Oo) + 62(^1 - ^0)' + o((^i - ^0)'), 
for 9 1 — 9q where bi = — sin 6*0, 
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and 

6^ 



COS02 = 1 - y + o{dl), for 02 0. 



Analogously to part i), we substitute these formulas into (16) to 
obtain 



2 



4 4 

32\ I „///) n \2\ 



= o{e'2) + om-e,y) + oiix-Xoy), 

for {9, A) -> {00, 0, Ao) where A = 6ai and B = > 0. Since ^'(A) ^ 
if -D(A) 7^ ±2 according to Lemma |4| D'(Xo) ^ and so A = Qai = 
3D'{Xo) 7^ 0. Thus again two spectral bands which touch at the point 
(^0, 0, Ao) have conical form near that point. The Bloch variety of 
operator H has conical singularity at {6q, 0, Aq). 
The same argument applies to (6*, A) = {—Oq, 0, Aq). 



Figure 6. Bloch variety of operator H in the free case 

2. When we turn on a small potential, since -D'(Ao) 7^ for Aq = 
D~^(±2/3), we can repeat the calculation in part iii), and find conical 
singularities at = {±60, 0) or (±6*0, ±7r) (at these points |-D(A)| = 
2/3). 

As it was shown in 22 , for almost every C°° periodic potential go 
on M, all the gaps Hill operator H^'^^ open (at the edges of the gaps 
D{X) = ±2). Since the spectrum a{H) has gaps iff a{HP^^) has gaps, 
this implies that all conical singularities in the free case with D{X) = ±2 
will generically split into two smooth branches and open a gap when 
we perturb the potential a little bit. □ 
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3. Proof of Proposition fTol 

Proof. 1. We will use the following notations a := cos 6*1 + 1, b := cos 6*2, 
then a e [0,2], 6 e [—1, 1]. Equation (|8| becomes: 

9x^ - X = a{3x + b). (18) 

Denote as the graph of function y = a{3x + b) and /i,/2,/3 
as parts of the graph of function y = 9x^ — x restricted to [—1,-1/3], 
[—1/3, 1/3] and [1/3, 1] correspondingly, / := /1U/2U/3. Then Ji,/2,/3 
are all connected. 

One can notice that the line /o,i intersects with Ii, 12,13 when x = 
— 1/3,0, 1/3 correspondingly while the line ^2,1 intersects with Ii, 12,13 
when X = —2/3,1/3,1 respectively. When slope 3a of the line la,i 
changes from to 6, the line 1 rotates from the line /o,i to the line /2,i 
around point (1,0). Thus, the line la,i intersects with each of Ji,/2,/3 
for all values a G [0, 2]. Applying the same argument for the line la-i, 
we also have that the line la-i intersects with each of /i,/2,/3 for all 
ae[0,2]. 

Besides, for a G [0,2], 6 G (—1,1), three lines la,i,la,b and la-i are 
parallel and la^ lies between the other two. Both lines la^i and la,~i 




Figure 7. The line li^i lies between /o,i and ^2,1 and thus 
has nonempty intersections with Ji, I2 and I3 



18 



NGOC T. DO AND PETER KUCHMENT 



intersect with each of Ji, /2, Is, and so the hne la b also intersects with 



each of Ji, I2, h- This means equation (18) has three real roots for all 



a G [0,2], b G [—1,1]. Thus all roots of equation ([8]) are real and so 
function F{6) has only real values for all ^ G 5. 




-1 -0.8 -0.6 -0.4 -0.2 



Figure 8. The line h/2,1/2 is parallel to h/2,1 and ^1/2,-1, 
thus intersects with Ji, /2, 13 



2. For each a G [0,2], b G [—1,1], the line la,b intersects with 
each of Ii, 12,13. By our notation, /i,J2,/3 are parts of the graph 
of function y = 9x^ — x restricted to [-1,-1/3], [—1/3,1/3] and 
[1/3,1] correspondingly. Thus, the ranges of functions Fi,F2 and F3 
are [—1, —1/3], [—1/3, 1/3] and [1/3,1] respectively. 

3. In what follows, we will find out when function F3 attains its max- 
imal and minimal values. The similar argument applies for functions 
Fi and F2. 

From part 2 we know that maximal value of F3 is 1 and its minimal 
value is 1/3. 



Plugging X = 1 into equation (18) we have a(3 + 6) = 8, which 
occurs only when a = 2,6 = 1, i.e. 9 = (0,0). So function Fi attains 
its maximum at (0, 0). 



Similarly, we plug x = 1/3 into the equation (18) to obtain a{l + b) 



0, i.e. either a = or 6 



-1. Now for each case when a = or 6 
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Figure 9. Graph of root function for equation dS 



solve equation (18), we will see that the biggest root of equation (18) 
is equal to 1/3 when a = 0, 6 G [—1,1] or a G [0,2/3],fe = —1. This 
means function attains its minimum at (±7r, 62) for 62 G [— vr, tt] or 
(ei,±7r)for 9^e[-n,-eo]U[9o,7i]. 

4. Let us denote the linear level set of function F{6) as L (if such a 
set exists). 

L := {{61,62) G B\p'l6i+p% = 2ko7T},po = (rf,p^) G Z^\{{0 , 0)} , ko G Z. 



For all 6 belonging to L, equation ([8]) has (at least) a constant solution, 
namely c. Then 

9c3-c= (cos^i + l)(3c + cose2), for all (^1,^2) G i^. (19) 
If P2 = 0, then the linear level set L = {{2kQiT/p^, 62), 62 G [— tt, tt]}. 



Since all values in (19) are constant except 62 changing from — vr to vr, 
the expression (19) is true only if cos6'i + 1 = 0. This would mean 
6i = ±71, i.e. L = {{±71, 62), 62 G [-vr,7r]}. 

In case P2 7^ 0, the linear level set L can be rewritten as 

P2 
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and so ( 19 ) becomes 



9c3 - c = (cos 01 + 1) (sc + cos ^^""^ ^ ' all ^1 ^ [-^' 

Since c is a constant and 6i runs from — tt to tt, we have 

9c3 - c = (cos 7r + 1) j^3c + cos -^—q J = 0- 

Thus by solving the equation 9c^ — c = 0, we conclude that constant c 



can be 0, 1/3 or —1/3. Plugging each value of c into (19), we can get 
all the linear level sets of F{6) as below: 

{{±n,92),92e[-n,n]}, 
m,±7r),ei G [-7r,7r]}, 
{(0i,±7r/2),0iG[-7r,7r]}, 
{(01,0), 01 G [-7r,7r]}. 

As a consequence, function -F(0) does not have any flat branches. □ 

4. Final remarks and acknowledgments 

(1) The lithographite structure G is not completely fiat (due to the 
presence of four bonds converging at some vertices) jsj. This, 
however, does not change the quantum network model that we 
study. Additionally, this structure is the least stable of the 14 
configurations studied in f^. It is, however, the easiest to study 
among graphynes. Indeed, complexity of the analysis grows 
with the number of atoms contained in a fundamental domain. 
This malces graphene the simplest (with just two atoms in an 
appropriately chosen fundamental domain) and the structure G 
of this worlc the next simplest, with three atoms. 

(2) The structure G has much less than a honeycomb symmetry, 
which does not prevent it from displaying Dirac cones. This 
happens also in various other graphyne structures (e.g. [jj). 



(3) As Figure 10 shows, the Dirac cones in this structure are highly 
anisotropic, indicating a very directional conductance. This di- 
rectionality effect, very sharply presented in the structure under 
consideration, has been noticed for several graphyne structures 
(e.g., [iSj)- This is one of the features malting graphynes fasci- 
nating. 



(4) 



In contrast to the density functional calculations of ff] for the 
structure G of this work, our results show no band overlap. 
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Figure 10. The local view of one of the stable Dirac 
cones. One notices that conductivity will be suppressed 
in one direction. 



Even though we and the authors of study different approx- 
imate models, they have the same geometry and thus are ex- 
pected to show coherent fetures. 

It is interesting to observe that computations for some other 
graphyne structures (e.g., # 10) in [7] do correspond well to 
such effects arising in the quantum network models. 
(5) The bound states arising in quantum graph models of graphene 
and graphyne (which cannot arise for the full dimensional peri- 
odic Schrodinger equations, e.g. 



I5 J2I]), should probably still 



suggest existence of some rather flat bands/ strong resonances 
(compare with the photonic crystal situation in [I6]). This is 
confirmed by computations in ff], except for the case of the 
lithographite structure G, which shows such flat bands in the 
quantum graph model, while these are absent in the results 
of(7|. 

(6) Due to the small symmetry group of the structure, the class 
of all invariant potentials is somewhat wider than the one we 
considered. Namely, it comes from two potentials: go on [0, 1] 
and qi (on [— 1, 1]Q It is somewhat harder technically to study 
this more general class of potentials, but the authors plan to 
address this issue in a future article. 

(7) The Corollary |8] and its analog for graphene (see [TT]) sug- 
gest that any non-trivial "obstacle" (potential) along the edges 



In this work, the potential qi is just the concatenation of two copies of Qq. 
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opens spectral gaps. For instance, the "supergraphene" struc- 
ture # 10 in |7], which differs from the standard graphene by 
presence of extra two atoms along each edge, is expected to 
and indeed does show gap opening (as well as very fiat bands 
at some gap edges, similar to the ones in our results). 
Besides studying more general invariant potentials, the aim of 
a future work is to consider spectra of nanotubes folded from 
the graphyne structure of this work. This, in particular, is the 



reason of the presence of the statement (4) of Proposition 10 
which will play significant role there. 

The authors express their gratitude to the referees for their substantial 
remarks and additional references. 
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